Generating coherent state of entangled spins 
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A coherent state of many spins contains quantum entanglement which increases with a decrease in the col- 
lective spin value. We present a scheme to engineer this class of pure state based on incoherent spin pumping 
with a few collective raising/lowering operators. In a pumping scenario aimed for maximum entanglement, the 
steady-state of N pumped spin qubits realizes the ideal resource for the 1 — > -y quantum telecloning. We show 
how the scheme can be implemented in a realistic system of atomic spin qubits in optical lattice. Error analysis 
show that high fidelity state engineering is possible for N ~ 0(100) spins in the presence of decoherence. The 
scheme can also prepare a resource state for the secret sharing protocol and for the construction of large scale 
Affleck-Kennedy-Lieb-Tasaki (AKLT) state. 

PACS numbers: 03.67.Bg,42.50.Dv,37.10.Jk,37.30.+i 



I. INTRODUCTION 

Coherent state in quantum mechanics usually refers to a 
specific type of quantum states with minimum uncertainty. It 
was first discovered in the context of oscillator field and has 
found wide applications in quantum optics O-0]. A quan- 
tum harmonic oscillator in coherent state most closely resem- 
bles the behavior of a classical oscillator. The notion was 
later generalized to spin systems JD-Hl]. In an ensemble of 
A*" spin-/ particles, the term spin coherent state (or atomic 
coherent state) is used to denote the states where the collec- 
tive spin J = ^2 n \ n has the minimum uncertainty JH]. Such 
states can be easily identified in the basis | J, fj,, A) which are 
eigenstates of J 2 and J 2 with eigenvalues J(J + 1) and fj, 
respectively, and A denotes additional quantum numbers to 
provide a complete set of labels. The collective spin J gets 
every value from NI down to (or i) in integer steps, and 
for each collective spin value J, the magnetic quantum num- 
ber /i = — J, — J + 1 . . . , J. One can easily show that the 
minimum uncertainty relation ((A J X ) 2 )((A J y ) 2 ) = \(J Z ) 2 
is satisfied for all extremal states \ J lt i = —J,X) in this ba- 
sis. Hence these states and their rotations generated by J are 
the spin coherent states (SCS) Hfl. Interestingly, two contrary 
characters coexist on these states: the most classical collec- 
tive spin behavior; and the fundamentally non-classical phe- 
nomenon of quantum entanglement. For every J < NI, there 
is a degenerate set of | J, — J, A) with identical collective prop- 
erties and distinct entanglements where the number of unen- 
tangled spins is upper bounded by j OJ. 

Preparation of SCS has been possible only in limited cases. 
The J = NI SCS, nondegenerate and unentangled, is ob- 
tained when all spins are fully polarized. Most experimental 
studies of spin squeezing start on this state. Schemes were 
also proposed to populate mixed state of singlets ( J = SCS) 
by collective pumping ifioll . and to select out singlet by pro- 
jective measurement in a scattering model Hill . Engineering 
a pure-state SCS of an arbitrary J value is a challenge but of 
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multi-fold significance. It is the sufficient condition for initial- 
ization into a decoherence free subsystem for robust quantum 
computation under strong collective decoherence 1121 . SCS 
of J <C NI are resources of large scale entan glem ent with 
potential uses in one-way quantum computation Il3lfl4ll . The 
ability to access a SCS of entangled spins also opens up a 
new realm for the study of the interrelation between collec- 
tive spi n behaviors and quantum correlations in a spin ensem- 
ble 

In this paper, we propose control schemes for engineering 
pure-state SCS of an arbitrarily specified collective spin value 
in a general spin ensemble. The schemes are based on inco- 
herent spin pumping of the N target spins and a set of an- 
cilla spins by a few (e.g. three) collective raising/lowering 
operators. The desired pure state is obtained with an N- 
independent probability by a single projective measurement 
on the steady state of the pumping, and the success rate ap- 
proaches 100% with O(10) cycles of pump plus measure. In 
a simplified pumping scenario aimed for maximum entan- 
glement, the steady state of N pumped spin qubits (with- 
out measurement projection) realizes the ideal resource for 
1 — > y optimal quantum telecloning lfl7ll . We show how the 
scheme can be implemented in the realistic system of atomic 
spin qubits trapped in optical lattice, where the collective spin 
pumping is realized by Stokes or anti-Stokes light scattering. 
Error analyses show that high fidelity state engineering is pos- 
sible for up to N ~ O(100) atomic spin qubits in the presence 
of control errors and decoherence. This is a concrete example 
of using simple and robust irreversible dynamics to prepare a 
desired complicated quantum state lfl8M23ll . The scheme can 
also prepare resource state for the secret sharing protocol l24ll . 
and for efficient construction of large scale AKLT state with 
applications in one-way quantum computation I25l - l27ll . 



II. GENERAL SCHEME 

Key to the state engineering approach by irreversible dy- 
namics is to design the dissipative controls under which the 
system saturates to the desired state vectors. We utilize here 
the spin pumping process which drives a spin system towards 
a mixture of all singlets connected to the initial state by the 
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pumping operators 111 OB . If a target spin ensemble is in sin- 
glet with a spin- J ancilla, its collective spin value must also 
equal to J. With a proper constraint from conserved quantum 
numbers, the singlet can be unique from which the desired 
pure-state SCS of the target spins can be obtained. 

The target spins are divided into two subgroups with col- 
lective spin ]a and js respectively, and the collective spin of 
2 J spin-i ancillas is denoted by j | g. When an inhomogeneous 

collective operator of the form J. = ca3a + cbJ% + 
acts on a SCS, the final state can be generally written as 



J+\Jt,-Jt,X) 



A, A' 



■A,-J T +1,A'). (1) 



where the first two quantum numbers in the kets denote the 
total spin and the z-component of Jt = + js + re- 
spectively. Calculation of the coefficients \ is straightforward 
by expanding the collective spin states in terms of common 
eigenstates of j\, j\, j B , j B , j 2 an d j% . We find that only the 
A = 0, ±1 transitions are allowed ill Oil , and the ratio between 
transition rates A = ±1 is 



./t + 1,A 



= (J T + 1)(2J T + 1) 



Jt + 1,A 
Jt,X 



(2) 



Consider the incoherent strong pump by J^, which results in 
a mixture of SCS of Jt, and the weak pump by the inhomo- 
geneous operator jf which then causes transitions between 
these SCS with the effective rate oc |x| 2 and the selection rule 
A = 0, ±1. From Eq. ||2), we can see the A = — 1 transi- 
tion is much faster than the A = 1 one between any such pair 
of states. Thus, the pump will saturate the target and ancilla 
spins to singlets of Jt where Jt is minimized. With the tar- 
get and ancilla spins initialized on the fully polarized state, 
the quantum numbers ja = NaI, jb = NbI and jp = J are 
all conserved by the pump operators. Only one singlet exists 
under this constraint: 
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where | J, fi, JAtJe) ab denotes eigenstates of (j^ 
3 a + 3 b w i tn eigenvalues J (J + 1) and p, respectively. 

Fig- Q2a) presents a simulation of the spin pump using 
the Lindblad master equation p = — | 5Zm=o(-^m-^mP + 
pL^Lrn - 2L m pL] n ), where L = y/AiJ^ and L m = 
y/AiJ+ for m = 1,2. Here we have chosen the inhomo- 
geneous raising operators 
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while other choices of coefficients ca,b,/3 lead to similar re- 
sults. For the simulated example, we set the spin pump rates 
Ah/Ai — 5000, and ja = 3b = jp = 5. After a pump 
time t p — 0.2A" 1 , \Sabp) is occupied by a population 



P(0) 
quire 



20%. For general values of ja, 3b and jp, we re- 
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FIG. 1: (a) Simulation of spin pump and repump by the collec- 
tive operators , and (see text). The target and the ancilla 
spins are in the fully polarized state at t = 0. Solid curve: pop- 
ulation on the singlet \Sabp)- Dotted curves: populations in the 
subspace of J = 1 (red), J = 2 (blue), and J = 3 (black) re- 
spectively. We assume the singlet is projected out at t — 0.4A~ , 
and hence the curve after is the repump dynamics, (b) Simulation 
of the simplified scheme for engineering the telecloning resource for 
N = 40 qubits. The populations on the singlet |0, 0, -j, -j), the 



spin-1 SCS |1,-1, 



-), and the spin-2 SCS |2, -2, 



-) are 



shown as functions of time using various pumping rates. 



which ensures the lowering operator to be applied much 
more frequently than the raising operator J. . The largest pos- 
sible value of ( Jr jf) is ~ [j a+ j b + j p) 2 while the smallest 
possible value of (J^Jt) is ~ 1. Thus Ah/Ai ^> (ja + 3b + 
jp) 2 is sufficient to ensure the condition in The steady- 
state population on \Sab@) is given by g(fc)] 1 = 20% 
where g(k) = (2k + (2t 2 + 3i + l)' 1 |l()|. The 

timescale to reach steady-state is ~ (ja + ]b + jp)" 1 N^ ■ 
The singlet can be selected out by projective measurement of 
J T or J T . If the measurement outcome is not singlet, the 
spins can be repumped to the steady state in a much shorter 
timescale ~ (j A + jB + jp^Ar 1 [Fig.QIa)]. The proba- 
bility of NOT obtaining \S A Bp) is reduced to 0.1% after 30 
cycles of measure plus repump. 

From the singlet \Sabp), further pumping by the target 
spin operator + j B bring the target spins to the desired 
SCS \ J,fi= -J,3a = NaI, 3b = N B I). Here the collective 
spin value J of the target spins is controlled by the number of 
the ancilla spins involved. Different choices of Na and Nb 
realize distinct pure SCS of the same collective spin value, 
which are fully symmetric under permutation of spins within 
subgroup A (or B). A and B can also be initialized with any 
jA < NaI and j B < NbI by applying the scheme first to the 
subgroups. Concatenation of the scheme can therefore realize 
pure SCS with more general permutation symmetries. 

SCS of the smallest collective spin values are most desir- 
able as a resource of entanglement. We consider the J = 
scenario of the above scheme (i.e. no ancilla spins) which uses 
two pump operators: the homogeneous j^ + j B and the in- 
homogeneous j\ ~ jg 



where A and B each contain A ^ 
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get spins. For A^/A; ^> N 2 , we find the steady state of 
the pumping p = P(J)\J, — J, f 7, § 7) (J, — J, f 7, f I| 
where P(J) = (2J 2 + 3J + 1)P(J + 1). This steady 
state is reached with a pump time t p w 3 + 1 ^ r ^ r/ A,^ 1 by our 
numerical estimation, and is largely a mixture of the sin- 
glet |0,0, f 7, f 7), the spin-1 SCS |1, -1, f I, f 7), and the 
spin-2 SCS |2, — 2, ^-7, -y 7), with some residue population of 
0.5% on the spin-3 SCS |3, -3, f 7, f 7). These states can be 
distinguished in a non-demolition way by measuring j z A +jg. 
A single cycle of pump plus measure thus ends up with one 
of these pure states which all have large scale entanglement. 
Fig- 02b) shows simulation of this spin pumping for a cluster 
of 40 spin qubits. 

The singlet |0, 0, ^ , ^) of N qubits turns out to be the ideal 
resource for universal optimal quantum telecloning ifrTll . If 
Alice holds subgroup A and each of her ^ associates holds a 
qubit in subgroup B, Alice can transmit identical copies of her 
unknown state cos ||0) + sin |e^|l) with a fidelity of Fo = 
to the -y associates using local operations and classi- 
cal communications (LOCC) I17I1 . By a single cycle of pump 
plus measure, the success rate to obtain this state is ~ 46%, 
which is a substantial improvement over the existing scheme 
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where the success rate is 111 111 . Most remarkably, all al- 
ternative outcomes by our scheme, i.e. | J, — J, ^) with a 
finite but small J, can also be used as quantum telecloning re- 
source under the same LOCC. Following the same procedure 
of Ref. Ql but replacing |0, 0, ^, f ) with | J, — J, f , f ), 
we obtain the telecloning fidelity which is then a function of 9 
[Fig. 12 a)], and it reaches the maximum value on the equator 
of the Bloch sphere. 

_ (3J + 4)A^ 2 + 4( J + 1)JV - 4J( J + 1)(J + 2) 
" 2{2J + 3)N 2 

Since Fj 1 ®* > Fo, better telecloning fidelity can be achieved 
with these finite J SCS in the presence of partial information 
(i.e. the range of the 9 value). For N 3> J, the telecloning 
fidelity averaged over the entire Bloch sphere approaches Fo 
[Fig. Ob)]. Thus, the mixed steady state of the spin pumping 
can be used as an equally efficient telecloning resource as the 
ideal singlet. 

A major cause of error for the state engineering is local 
spin decoherence process. If each spin loses its phase co- 
herence with a rate 7, a total leakage of ~ Njt p ~ 7/ A, 
out of the desired subspace is accumulated in the entire du- 
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of the state preparation. High fidelity 
state engineering thus requires: 7 <C Aj. This is confirmed 
by numerical simulation for a cluster of N = 8 spin qubits 
where we have added pure dephasing processes described by 
Lindblad operators y/2^yl* for all spins [Fig.0a-d)]. We also 
studied the effect of errors from system parameters. For the 
simulation presented in Fig.|3e-h), spins are pumped instead 
by + E B and ~\ - S+ where E ± = £ n (l + r?„)7±, 
rj n being a random error between rj and —77. The figure of 
merit is reasonably good when the error amplitude rj < 10%. 
Moreover, by pumping with the operators U(j A + jg)U' 
and U(j\ - Jb)U^ where U = 11™ ex P(*#n-^n)> the state 
U I J, - J, f 7, f 7) is obtained instead of | J, -J, f 7, f 7>. 
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FIG. 2: (a) Fidelity of 1 -> 20 telecloning of a state cos §|0) + 
sin|e l< ^|l) using the above singlet (solid), spin-1 SCS (dashed), 



spin-2 SCS (dotted), and the spin-3 SCS |3, 



-) (dot-dashed) 



as the resource respectively, (b) The fidelity of 1 — ¥ telecloning 
(averaged over the Bloch sphere) using the above resource states re- 
spectively. 



Namely systematic phase errors in the collective pumping op- 
erators do not affect the entanglement, and single spin ro- 
tations about the z-axis can be deliberately encoded in the 
pumping. 

When A and B each contains 2 qubits, the resultant sin- 
glet |0, 0,1,1) by our scheme is the 4-qubit AKLT state 
P2s\S)i2\S)s4, where \S)ij stands for the singlet of qubit 
i and j, and Pjk is the projection operator to the triplet 
subspace for qubit j and fc j25ll . Its optical an alog has 
been used to demonstrate four-party secret sharing []24[, and 
measurement based single qubit rotation 1I26I1 . This state is 
also an efficient element to construct a large scale AKLT 
state as schematically illustrated in Fig. |4] Consider two 4- 
qubit clusters in P 2 3\S)i2\S) 34 and P67\S) 56 \S) 7S respec- 
tively, by measuring the parity of atom pair {4, 5}, the spin 
configuration of this pair will be projected to either the sin- 
glet or the triplet subspace J28|. With 75% probability, the 
measurement outcome is triplet and an 8-qubit AKLT chain 
^23^45^6715)12^)341^)561^)78 is obtained. The rest 25% 
probability will give 7 , 23-P67|5)i2|5)3 6 |5)7 8 <8> |5) 45 where a 
6-qubit AKLT state is obtained. With our scheme as an effi- 
cient source of 4-qubit AKLT states, a long AKLT chain can 
thus be constructed. 



III. APPLICATION TO ATOMS IN OPTICAL LATTICE 

Here we apply our scheme to cold atoms of a typical A-level 
structure which are trapped in optical lattice. This system has 
been widely explored in various schemes for quantum infor- 
mation processing. The two lower energy levels are used to 
represent the spin qubit (Fig. |3a)). The atoms trapped in an 
optical lattice can be first initialized to the | \) state and loaded 
into a Fabry-Perot cavity (Fig. |5jb)). For simplicity, we as- 
sume the optical lattice constants in both directions equal to 
the wavelength of the cavity mode and all atoms are at the 



4 




0.05 0.10 



FIG. 3: Merits of state engineering in presence of spin decoherence (a-d) and system parameter errors (e-h) for engineering the telecloning 
resource for N = 8 qubits. Au/Ai = 100. (a) and (e): Probability of obtaining singlet (blue triangle), spin-1 SCS (green square) and spin-2 
SCS (red diamond) at t p = 0.75A" 1 . (b) and (f): Bipartite entanglement E(A\B) in the obtained states measured with the logarithmic 
negativity, (c) and (g): Entanglement between one qubit in A and another in B. (d) and (h): Fidelity of the obtained states with the target states 
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• atoms — singlet 



projection to triplet 



FIG. 4: Increasing the length of a AKLT chain with the resource of 
4-qubit-AKLT-state. By measure the parity of the two atoms, the 
n-qubit AKLT chain either becomes a (n+4)-qubit (lower left) or a 
(n+2)-qubit AKLT chain (lower right). See text. 



peak of the cavity field. Driven by two lasers of frequency 
uj c ± u> z with u! c being the cavity resonance and oj z the spin 
splitting, the cavity-assisted Raman process lowers/raises the 
spin state. The Raman processes in the large detuning regime 
(g 2 , fi± -C A 2 ) can be described by the effective Hamiltonian 



h.c, 



(5) 



where Sl± is the Rabi frequency of the two pumping lasers, 
g the atom-cavity coupling and A± the detuning. a\ creates 



laser wavevector k. For the setup shown in Fig. |5jb) where 
the pump lasers are perpendicular to the cavity axis, 'blue' 
laser with 6\ = ^ realizes the homogeneous operator , 
and 'green' laser with cos 6% = | (or |) realizes the inhomo- 
geneous collective operator (or J"^~) where subgroup A, 
B and ancilla are represented by blue, red and green spheres 
respectively. 

Projective measurement for selecting out the singlet state 
can be realized in the same setup. Applying a 'blue' and a 
'green' laser with both # = § and comparable Rabi frequency 
fi_ ~ + realizes the homogeneous raising and lowering 
operators on the spin qubits. If the system is in finite J 
state, then Raman scatterings are allowed and we shall observe 
continuous cavity photon emission when jjt and pump 
the spins. When the system is in singlet, Raman scattering is 
forbidden since both j£ and annihilate the state and there 
will be no cavity photon emission. 

The Raman scattering rate by a single atom is A^u = 
PrSl^/A^ with r being the spontaneous emission rate of 
atomic excited state in vacuum and P the cavity induced 
enhancement factor (Purcell factor). Consider the Cs atom 
(j- = 2.6 MHz) in a typical Fabry-Perot cavity with mode 
volume of 10 4 /^m 3 and quality factor 1.7 x 10 7 , which cor- 
respond to P sa 80, ~ 45 MHz and cavity decay rate 
20 MHz fish. A ps 15 MHz can then be achieved 



-'k-rj jji realizes various 
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with 
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150 MHz and 



40 MHz [301. The col- 



a cavity photon and J (k) = 
collective spin raising/lowering operations by controlling the lective Raman scattering rates shall satisfy Aj(J{" 
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(a) 



atom j 




Ah( J£ j^) < k (Tfl] , the last inequality is to ensure the emis- 
sion of cavity photon is spontaneous. For N atoms in cavity, 
the matrix element (J^J^) ({JiJ±) ) is ~ N (~ AT) in the 
neighborhood of the polarized initial state and ~ 1 (~ 
in the neighborhood of the target singlet state. We can thus 
use a Ah ~ 15 MHz. Correspondingly, Aj shall be ramped 
down from an initial value of ~ 1 MHz to the steady state 



value 



N- 



MHz along with the spin pumping. Since atom 



in optical lattice can be of an ultra-slow spin decoherence rate 
■2- ~ 1 — 25 Hz Olll . the condition 7 <C A, can be satisfied 
for N ~ 0(100) qubits. 



FIG. 5: (a) The level structure of the atoms. Red arrow represents the 
cavity field, and blue and green arrows denote two laser fields, (b) 
Collective spin pumping of atoms realized by cavity assisted Raman 
process. Atoms are trapped in optical lattice and loaded into a Fabry- 
Perot cavity. The inhomogeneous coefficients are controlled by laser 
emission angles 9i and 62- 
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